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Abstract. - 

Inelastic grains in a flow under gravitation tend to collapse into states in which the relative 
normal velocities of two neighboring grains is zero. If the time scale for this gravitational 
collapse is shorter than inverse strain rates in the flow, we propose that this collapse will 
lead to the formation of "granular eddies" , large scale condensed structures of particles moving 
coherently with one another. The scale of these eddies is determined by the gradient of the strain 
rate. Applying these concepts to chute flow of granular media, (gravitationally driven flow down 
inclined planes) we predict the existence of a bulk flow region whose rheology is determined 
only by flow density. This theory yields the experimental "Pouliquen flow rule", correlating 
different chute flows; it also correctly accounts for the different flow regimes observed. 



Introduction. - Flows of hard granular systems are ubiquitous in nature and technology, 
yet are still poorly understood 0. Granular systems typically have a twofold separation of 
energy scales: the typical energy of a particle is determined by gravity or some other body force 
(in a few instances by initial conditions), and is much larger than the thermal scale fceT, yet 
much smaller than the scale required to appreciably deform the particle. Despite the smallness 
of ksT on the scale of granular energies, many treatments use a pseudo-temperature connected 
to the random part of the kinetic energy of a particle. Such treatments often link granular 
phenomena to the kinetic theory of gases. The "granular gas" has an intrinsic rheology, and 
is driven by the external forcing. 

One of the pioneering treatments of this rheology was by Bagnold, who discussed chute 
flows, the gravitationally driven flow of a granular material down an inclined surface 0. It 
is simplest to consider a flow of constant, fixed depth H, with the average velocity of the 
particles parallel to the free surface. The particles are spheres of monodispersed mass M and 
radius R. We choose axes such that the direction of flow is x, the direction perpendicular 
to the free surface of the flow is z, and the direction parallel to vorticity is y (see Figure |l|). 
The shear stress o~ xz in such a flow is communicated by particles at slightly differing depths, 
whose velocities differ if d z v x is non-zero. We expect the momentum transfer communicated 
by collisions between particles at different depths to be of the order of 
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Fig. 1 - (a) Chute flow is flow down a surface inclined at an angle 9. The a;-axis is chosen parallel 
to the flow, the z-axis perpendicular to the free surface. The y-axis is parallel to the vorticity of the 
flow, and is directed out of the page, (b) In the granular eddy picture, the motion of the particles is 
regarded as a superposition of the translation and rotation of granular eddies. 



Ap = MRd z v x . (1) 
Furthermore, these collisions will occur at typical intervals of the order of 



At=(d z v x )-\ 

from which we conclude that the typical collisional stress will be 
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In a steady state flow down a surface inclined at an angle 9, the xz shear component of 



the stress tensor is determined by gravity to be 



<Jxz = pgz sin ( 



(4) 



with p the (local) mass density, which we here assume to be independent of z (we will return 
to this point below.) We are measuring the depth of the pile z from the free surface, at which 
z = 0. 

If a xz oc z oc (d z v x ) 2 , then since v x = at the base z = H if the boundary pins the flow, 
we immediately obtain 



d z v x 



-At 



defining the coefficient Asag, or 



v x (z) 
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(5) 



(6) 



While there have been a variety of authoritative experimental studies of chute flow g , as 
well as intricate theoretical discussions of the rheologies to be expected on general grounds B , 
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we have been particularly inspired by the recent work of Pouliquen . Pouliquen studied the 
behavior of chute flows as a function of inclination angle 9 and height of the flow H . He found 
that for small values of 9 or height H , no flow took place. 

With the increase of either 6 or H such that an angle of repose line 9n(H) was passed, 
a region of steady-state flow was entered. Finally, for values of 9 above a maximum 9m, the 
flows continuously accelerated, and no steady-state flow was observed. 

The dominant observational fact about the steady-state flows is the "Pouliquen flow rule," 
which connects the vertically averaged velocity u of a flow of height H with the height H stop 
at which flow ceases for a chute of that inclination 9. [The angle of repose 9n(H) is the inverse 
of the function H stop {9).\ The Pouliquen flow rule gives a scaling form for the average velocity 
u of the flow, 



f 

of 



H 

v x (z)dz, (7) 



= 0-5?-. ( 8 ) 

' ' stop 



and accounts well for experimental data with (3 = 0.136. 

The scaling u oc H 3 ^ 2 in the Pouliquen flow rule is consistent with the Bagnold rheology. 
But the Pouliquen flow rule also connects the coefficient A^ ag (9) with the thickness of the pile 
at that inclination below which flow arrests, which would not be expected from the Bagnold 
point of view. 

Note that we would have obtained dimensionally the Bagnold result for the rheology had 
we claimed that the stress should obey 

<7xz = lid z v x ~ pR 2 (d z v x ) 2 , (9) 

where we have made the substitution for the viscosity /i ~ pR 2 d z v x on grounds that a granular 
flow has no other obvious local length or time scales than R for the length scale and (9zU K ) _1 
for the time scale. The heart of the Bagnold approach thus lies in the assumption that these are 
in fact the only local scales. Note that the gravitational constant g does not figure directly 
in either of these scales. However, the Pouliquen flow rule implies that this rheology docs 
depend both upon g and upon the thickness of the arresting pile H stop , which is hardly local 
information. Thus the Pouliquen flow rule appears to be inconsistent with any assumption of a 
purely local rheology comparable to that of a granular gas This, and other considerations, 
have motivated some authors to build non-local models for the rheology H . 

The broad features of Pouliquen's conclusions have been confirmed by a series of numerical 
studies in which these authors have participated Jlffl. For relatively thin piles, the Bagnold 
rheology breaks down (as also seen in experiment, [g]), but the thicker piles show a Bagnold 
rheology and obey the Pouliquen flow rule, albeit with a slightly larger value of (3 (The 
crossover is examined numerically in ||). However, the assumption of the Bagnold or granular 
gas approach that the stress is mostly transmitted through collisions seems not to be true in 
these numerical studies; stress seems instead to be transmitted primarily through relatively 
long-lived contacts between particles. The density in the interior of the piles is independent 
of depth, consistent with the assertion made in Eq. (^) . 

In this treatment, we eschew granular gas approaches, and we do not assume the exis- 
tence of any rheology independent of the gravitational character of the flow. We show that 
gravitation combined with particle inelasticity is able to dissipate a significant fraction of the 
system's kinetic energy over time scales short compared to the inverse strain rate. Given this 
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fact, it is natural to posit the existence of "granular eddies," gravitationally collapsed networks 
of particles, which move coherently and whose properties determine the flow rheology. The 
rheology that follows from this picture agrees with the Pouliquen flow rule, and also gives a 
simple explanation for the different flow regimes observed by Pouliquen. 

In this report, we first analyze the phenomenon of gravitational collapse for inelastic parti- 
cles. We then introduce the granular eddy picture, and relate the eddy size to flow properties. 
We specialize to the case of chute flows, determining the rheology, and accounting for the 
principal features of the observed phenomenology of these flows. 

Gravitational collapse. - Consider an inelastic ball with a coefficient of restitution of e, 
bouncing on a rigid horizontal surface. It is elementary to show that if its normal velocity at 
first impact is vq, then after a finite time r gc it will come to rest, with 

T gc = "j ■ (10) 

g 1-e 

(A similar result obtains if we take a more realistic ball with a Hertzian contact force and a 
visco-elastic dissipation; here we restrict ourselves to the simplest case.) 

Now consider a particle in a granular flow. Suppose that r gc is short compared to the 
time scales for its neighboring particles to rearrange themselves. This latter time scale is the 
inverse of the local shear rate 7; for chute flow 7 = d z v x . In this short "collapse time" regime, 
an individual particle will rapidly come to rest upon the particle or particles beneath it in the 
flow; it may then roll about, possibly arriving in a stable configuration with three particles 
beneath it. However, its motion is strongly constrained by its gravitational collapse, and its 
motion is strongly correlated with that of the particles with which it is in contact. 

Thus we can envision large aggregations of particles coming into existence, each of whose 
motions with respect to its neighbors is at most of a rolling kind. Influenced by pictures of 
turbulence as a superposition of eddies, we term these aggregations "granular eddies" (see 
Fig.0). 

Granular eddy size. - Consider an eddy of scale £, whose center of mass is at a position 
Zq. The local average velocity can be expanded as 

v x (z) = v x (z ) + (z- z Q )d z v x (z Q ) + i(z - z ) 2 d zz v x (z ) H (11) 

While the first and second terms can be matched by an eddy whose center of mass moves 
at a velocity v x (zq) and which rotates at an angular velocity to = d z v x (zo), the third term in 
this series is incompatible with the rigid-body rotation of an eddy. We write this incompatible 
velocity at the eddy boundary as 

v ic = \e dzzVx . (12) 

Supposing this incompatibility is relieved by internal strains of the eddy on the scale t of 
the eddy itself, the time scale corresponding to variation of these "incompatibility strains" is 

I (i \~ x 
r ic = — = -d zz v x . (13) 

Vic V 2 / 

On this time scale the environment of particles at the boundary of an eddy will inevitably 
change as that eddy conforms with the surrounding flow. 

Since the radius £ determines the location of the outer perimeter of the eddy, then we 
anticipate that particles associated with neighboring eddies will typically be colliding with 
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the eddy with relative velocities comparable to v c — 2lo£. The characteristic gravitational 
collapse time associated with these collisions will be 



t,„ = (14) 



9 1-e 

so that the criterion for the particles at the eddy surface to be able to gravitationally collapse 
before their environment is altered by incompatibility strains is 



' gc 



ud z v x 



t ic \ g 1 - ej \2 



d zz v x < a, (15) 



where a is an unknown numerical constant of 0(1). The maximum value of £ consistent with 
this relation is determined by 

t\ [d z (d z v x ) 2 ] = ~ag X —^-. (16) 

We expect that this maximum value will set the scale of the eddies, since eddies smaller 
than this size will tend to grow as more and more particles collapse onto their perimeters, and 
eddies larger than this sizewill lose particles from their perimeters. 

Phenomenology of chute flow. - Let us use dimensional analysis to define an effective 
"viscosity length scale" £ v by inverting the Bagnold scaling relation given by Eq.(||) after 
substituting this new length scale instead of R. This yields 

a xz = p£l(d z v x ) 2 . (17) 
For chute flow, this can be combined with Eq.(||) to give 



' ■ - \l & 

We now make a different scaling assumption than that of Bagnold, which is that the length 
scale appearing in Bagnold's argument should be set by the eddy size £ e instead of the particle 
size R, i.e., 

£ = i»(l + &~ + ---), (19) 

where the unknown numerical constant b accounts for the leading order finite-size corrections 
due to the existence of the particle size R. In other words, we are assuming that (with the 
exception of these corrections in R), this length scale is the unique length scale determining 
the bulk rheology of the granular dispersion. Then £ e and d z v x are jointly determined by 



simultaneous solution of Eqs. (|16|),(|18|) and (19) 



If the angle of inclination is too small, there is no solution, in particular, for 



< 6 R = sin" 1 ( h- — - ) . (20) 



For 9 > Or, we have 
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which fixes the eddy size £ e . Finally, we find 

d z v x = -A Bag yfz (22) 

with 



sin6 R y/gsm6 R ( smO \ 

^Bag = 7 = = - — 7 1 J . 

te bR \sm6 R J 

We have up to now ignored the question of the density of the flow. Let us suppose that the 
eddies themselves have a fixed density po, independent of their size. Then the medium as a 
whole can have a density that differs from this ( "free volume" ) only due to a presumably lower 
density in the regions, of typical scale R, that separate different eddies from one another. We 
conclude that 

P = Po (l ~ + ■ ■ •) (24) 

or 

e^ = U^f-i), (25) 

po b \sm6 R J 1 

where c is yet another unknown numerical constant. Note that for chute flow, the eddy size 
£ e given by Eq.([2l|)is independent of depth (z), consistent with the assumption of a depth- 
independent density p. 

Of course, the eddy size £ e is not entirely unrestricted. When the eddies get too large to 
be accomodated within the height H of the flow , i.e., for £ e ~ H/2, we expect flow arrest 
to occur since all particles in the flow become connected to the bottom surface through the 
contact network. This gives the thinnest flowing pile at a given angle 8 in terms of a new 
constant d ~ 1 , 

sin 9 ' 1 



H atop {9) = d£ e {9) = bdR —- -I) , (26) 



smt 



or, equivalently, the lowest possible angle of stable flow at a given pile height H, 

9 n {H) = sin" 1 (^sin 9 R + bd^j ^ . (27) 

On the other hand, gravitational collapse ceases to stabilize the flow when £ e ~ R, corre- 
sponding to an upper limit of stability 

M =sin- 1 (esin(0 iI )), (28) 

which is independent of the flow height H. Here e ~ 1 is a further unknown proportionality 
constant. Note that we can now re- interpret 9 R as the limiting value of the angle of repose as 
the thickness of the flow H — * oo. 

Since this is a scaling theory, it is not possible to make a quantitative comparison between 
these predictions and numerical results such as those of Ref. flC] . However, for 6 — 6 R <C 1, 
the dependences on tilt angle for a thick flow (H ^> R) such as d z v x — — A^^-Jz with 
oc (9 — R ), and po — p cx (9 — 9 R ) are borne out by these numerical results. Also, the 
Pouliquen flow rule Eq. (W) is recovered with 
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H 




Fig. 2 - The granular eddy picture predicts that the angle of repose line separating the region of no 
flow from that of stable flow depends on H. The angle of repose line approaches a fixed angle 9 a 
as H — > oo. The maximum angle of stable flow 9 m does not depend on H in this picture, except at 
small values of H where the stable flow regime disappears. 



0=— VstiJ^. (29) 
5 

Finally, the form of the phase boundaries (see Fig. ||) is in agreement with both numerical 
and experimental results. 

In this Letter we have addressed the form of the bulk rheology for chute flows; our con- 
clusions regarding this rheology should hold in portions of the pile for which the computed 
eddy scale I is less than the distance to the boundaries. Clearly there will be both upper 
and lower boundary layers in which this is impossible. We have not addressed the structure 
of these boundary layers. Although chute flow does seem to have a bulk-rheology dominated 
regime, this may not be the case with all flow geometries. In some flows the structure of the 
boundary layers may dominate in determining the characteristics of the flow. 

* * * 

We are grateful to P.M. Chaikin, J. Landry, and L. Silbert for helpful discussions. 
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